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Abstract 

This paper presents a single frame algorithm for the spin-axis orientation- 
determination of spinning spacecraft that encounters no ambiguity problems, as well as a 
simple Kalman filer for continuously estimating the full attitude of a spinning spacecraft. The 
later algorithm is comprised of two low order decoupled Kalman filters; one estimates the 
spin axis orientation, and the other estimates the spin rate and the spin (phase) angle. The 
filters are ambiguity free and do not rely on the spacecraft dynamics. They were successfully 
tested using data obtained from one of the ST5 satellites. 


I. Introduction 

In the early days of space exploration, the use of spinning satellites was prevalent for 
spacecraft (SC) stabilization [Reference 1, Chapters 10, 11]. In that era only batch algorithms 
were used in order to determine the attitude of the spinning satellites. Starting in the late 
1970s, the focus has shifted from spinning satellites to three-axis stabilized SC [Reference 1, 
Chapter 12]. Considerable effort was invested in devising accurate algorithms for attitude 
determination (AD). In particular, a variety of recursive AD algorithms were introduced. As a 
result, spinning SC development and their resulting ground system development stagnated. In 
the 1990s, shrinking budgets made spinning SC an attractive option for science. The attitude 
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requirements for recent spinning SC are more stringent and the ground systems must be 
enhanced in order to provide the necessary attitude estimation accuracy, and yet suitable 
recursive algorithms for spinning SC did not exist. Therefore when the use of spinning SC re- 
emerged, efforts were made to develop such algorithms. Baker 2 , for example, developed a 
Kalman filter (KF) based on a dynamic model presented by Markley, Seidewitz and 
Nicholson 3 . The attitude was represented by Euler angles. The first derivatives of the states 
were nonlinear (trigonometric) functions of the states themselves. Simplifying assumptions 
had to be adopted in order to use the dynamics model in an extended KF. Sedlak 4 used 
Markley Variables 5 to describe the. SC attitude dynamics.. These variables are slowly varying 
which facilitates the filter state tracking and estimation, but the models which have to be used 
in the KF are quite complicated. Bar-Itzhack and Harman used a pseudo linear filter 6 to do the 
same 7 . The philosophy that governed the newly developed recursive filters for AD of spinning 
SC was an extension of the concepts on which three-axis stabilized AD algorithms were 
founded. Accordingly, other than [7], there was no separation between the spin axis 
orientation states and the spin angle states. Thus, the slowly changing dynamics of the spin 
axis orientation was combined into one dynamics model that included the fast changing spin 
(phase) angle. 

The present algorithm is based on the premise that the parameters which describe the 
direction of the spin axis orientation in inertial space vary slowly even when the SC nutates 
and precesses. The spin (phase) angle, on the other hand, changes fast but stays almost at a 
constant rate per a single revolution and is decoupled from the other axes. Iri fact, this is the 
classical approach to spin-axis orientation determination (ORD) [Reference 1, Chapters, 10, 
11]. (A good exposition of the difference in approach to three-axis AD and spin-axis ORD is 
presented in Reference 8.) This realization enables the decoupling of the recursive AD 
algorithm into two simple low-order filters that are independent of the SC dynamics. 
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There are two approaches to spin-axis ORD. One relies heavily on the solution of 
trigonometric functions [Reference 1, Chapters 10, 11]. The other approach is a vectored 
approach 9, 8 . This is the approach adopted in the present work; however, unlike in References 
8 and 9, here we develop two recursive algorithms, one for obtaining a single frame solution 
and the other is a novel Kalman filter for time varying ORD which is based on multiple 
measurements performed at different time points. Also, here the components of the spin axis 
are found as projections on the axes of the Geocentric Inertial Coordinates (GCI) rather than 
projections on the measured directions and their cross product, as presented in Eq. (1) of [9], 
Eq. (11 -3e) of [1], and Eq. (2) of [8]. In the present work we are not concerned with the 
measurement techniques. This topic can be found in other works [see e.g. References 1 and 
10 ]. 

As is well known, when only two vector measurements are available for spin axis 
ORD, there exist two possible solutions [see e.g. Reference 1, Chapters 10, 11 and Reference 
8]. The cause of this ambiguity is explained and a solution is proposed, which does not rely on 
cumbersome spherical geometry solutions. 

In the following section we discuss the geometry of the ORD problem. The ambiguity 
problem generated by the existence of two possible solutions is explained in Section III 
whereas in Section TV we present a simple vector solution to the ambiguity problem. Section 
V presents a simple low order Kalman filter (KF) for the spin axis ORD, and an even simpler 
one for estimating the spin (phase) angle. Results are presented in Section VI, a discussion of 
these results is given in Section VII, and the conclusions from this work are presented in the 
last section. 
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n. Connections between the Spin Axis Orientation and Vector Measurements 
Consider Fig. 1 where the sun sensor measurement is expressed by the components of 
the unit vector in the sun direction, resolved in the GCI coordinate system. These components 
are s lx , s Iy and s lz . Similarly, the normalized three-axis magnetometer (TAM) measurement 

is expressed by the three components, m lx , m [y and m lz (the TAM vector, m , is shown 
in the figure but its components are not shown). We want to find the direction of Z b in the 
GCI coordinates expressed by its components along the coordinate axes. When we know 
these components, we can certainly express the direction, of Z b by the angles a and p , if 

required. Denote the components of Z b in the inertial coordinates, GCI, by x, y and z. In the 
filter that will be presented later, we estimate Z b where Z b = [x y z] , 



Fig. 1: The Geometry of the Spin Axis and the Sun Vector 
Since the sun angle, cp s , is measured, we can calculate its cosine. Let us denote this 
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cosine by U s ; that is 


( 1 ) 


U s =cos<p s 

We note that the cosine of <p s , which is the angle between the sun vector and Z b , is nothing 
but s bz . On the other hand the cosine of cp s is equal to the dot product of Z b and s ; hence we 
can write 

U s = SfcX + s Iy y + s b z (2) 

Like with the sun sensor measurement, the cosine of the angle between the normalized TAM- 
vector and Z b is simply m bz (see Fig. 2); that is, 

* 

cos(p m =m bz (3) 

Denote this cosine by U m 

U m =cos(p m (4) 

Like with the sun sensor, we know that Z b • m = coscp m , hence 

U m= m txX + m Iy y + rn Iz z (5) 



Combining Eqs. (2) and (5) into one matrix equation yields 
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III. The Ambiguity Problem 

Equation (6) can be satisfied by two solutions. This is because, as is well known 9 ’ 1 ’ 8 , 
the spin axis can lie along two different lines. As shown in Fig. 3, these lines are formed by 
the intersection of two cones; namely, the sun cone and the magnetic field cone. These cones 
are described as follows. The main axis of the sun cone is the sun line* s , and the main axis of 
the magnetic field cone coincides with the normalized magnetic field vector, m . 



Fig. 3: The Geometry that Depicts the Two Possible Solutions. 



0 


Fig. 4: Upper View of the Two Possible Solutions. 
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The sun cone is generated by a rotation about the sun direction of a line that forms the cone 
half-angle (p s with the sun line. Similarly, the magnetic field cone is formed by a rotation 

about the magnetic field line of a line that forms the cone half-angle (p m with the magnetic 
field line. The two possible solutions are the line from 0 to p, and the line from 0 to p 2 (the 
lines themselves are not shown in the figure). An upper look at points p, , p 2 , s and m is 

presented in Fig. 4. The existence of two possible solutions can be demonstrated through the 
following example. 

Suppose that Z bI , the spin axis unit vector expressed in the GCI coordinates, is as follows 


Z w = 


0.612 

0.354 

0.707 


(13.a) 


Also suppose that 



'-0.176' 



' 0.162 ' 

s > = 

0.44 

(13-b) 

m, = 

-0.566 


0.88 _ 



0.808 _ 


where s, and m, are, respectively, a unit vector in the sun direction, and a unit vector in the 

direction of the magnetic field, both expressed in the GO coordinates. Then the matrix H , 
which is embedded in Eq. (6), and is defined as 


takes the value 


Six Sty s lz 


L m ix 

m i y 

™IzJ 

-0.176 

0.44 

0.88 

0.162 

-0.566 

0.808 


(13-d) 


(13.e) 


In this case U s = Z b • s, == 0.67 and U ffi = Z b • m, = 0.471 . The following two solutions 
Z M1 =[-0.743 -0.098 0.662] and Z bI2 = [0.612 0.354 0.707] satisfy Eq. (6). Indeed, 
Eq. (6) takes the following two correct forms 
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'-0.743 

0.67 


'-0.176 0.44 

0.88 1 

-0.098 

— 



0.471 


0.162 -0.566 

0.808J 

0.662 








'0.612' 

0.67 


-0.176 0.44 

0.88 ' 

0.354 

— 


_0.471_ 


0.162 -0.566 

0.808 

0.707 



L. 

* J 


(13-f) 


(».*) 


A way to resolve the ambiguity problem is presented next. 


IV. Ambiguity Resolution 

While other ways to solve the ambiguity problem also exist [see e.g. Ref. 8], we 
chose a rather simple approach to resolve the ambiguity problem, and indeed to solve for Z M , 
which avoids tedious spherical geometry calculations. Our solution to the problem is as 
follows. 


Let s t and nij be, respectively, unit vectors in the direction of the sun and the 
magnetic field resolved in the GCI system. Similarly let s b and m b be the same in the body 
system. Let the transformation matrix from the body to the GCI coordinate system be denoted 
by D, b , then obviously 


let 


then 


[ Sj | m, | Sj x m, ] = D; [ s b | m b | s' b xm b ] 

C = [ Sj | m, | s, xm,] (14.b) B = [ s b | m b | s b xm b ] (14.c) 


Dj b = CB 


-1 


(14.a) 


(14. d) 


Now Z bb , which is Z b , resolved in body coordinates, is given by Z bb T = [0 0 l] , and since 

(14.e) 


then 


2<bl — Dj z bb 


v=<*» 


(14.f) 
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where d b 3 is the third column of D b . Note that the two vector measurements have to be taken 
at certain time points; namely, at times when the sun sensor acquires the sun. 

The method proposed here is demonstrated through the following example. Let Z bI , 

Sj and m, be as before. Suppose that the SC is oriented such that 



' 0.227 " 




0.875" 

S b = 

-0.706 

0.67 

(15.a) 

and 

m b = 

0.113 
0.47 1_ 


(15 b) 


where s b is s resolved in the body coordinates and m b is the m vector also resolved in the 
body frame. Computing the C , B , and D b matrices defined', respectively, in Eqs. (14.b, c 
and d) yields 



"-0.176 

0.162 

0.854" 



" 0.227 

0.875 

-0.408" 

C = 

0.44 

-0.566 

0.285 

(15.c) 

B = 

-0.706 

0.113 

0.48 


_ 0.88 

0.808 

0.028 



0.67 

0.471 

0.644 _ 


-0.242 



-0.768 


0.593 


0.753 0.612 

-0.534 0.354 
-0.385 0.707 


(15. e) 


It is seen that the third column of D, b is Z b , given in Eq. (13. a) which is Z bI2 , the second of 
the two solutions, found before, of Eq. (6). It should be noted that D b can be found using the 
TRIAD algorithm 11,12 where unlike in Eq. (14.d), there is no heed to invert a matrix. We 
chose to use the present method for computing D b because TRIAD is a more elaborate 
routine whereas, being a 3x3 matrix, the inverse of B can be computed analytically. 


V. A Novel Kalman Filter for Determination of the Attitude 

The filter that is presented here consists of two linear reduced-order filters. The first 
filter estimates the spin axis orientation, Z bI , whereas the second filter estimates the spin 
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(phase) angle. Once the initial estimate of Z w is close to its correct value, the filter 
recursively produces a better estimate of it. The purpose of the discussion of the ambiguity 
and its resolution presented in the preceding sections was aimed at supplying the correct 
initial estimate, free of ambiguity. 

V.l A simple Low-Order KF for spin axis determination 
Recall Eq. (6) 


U s _ S I X S Jy S Iz 

uj- rn IX m Iy m Iz 


Since U s and U m are the result of measurements, we add to the last equation some 
appropriate zero-mean white-noise components v s and v m , and obtain the measurement 
equation 


U s _ S Ix S Iy S Iz + V s 

U m m ix m Iy m lz v n 

z 


(16.a) 


Since between measurements the direction of the spin axis of a spin stabilized SC does not 
change much, it is appropriate to model the dynamics of its components between 
measurements as a Markov process 13 . That is 


-1/t„ 0 6 Yx 1 fw. 


y = 0 -l/t y 0 y + w y 

z 0 0 -1 /t, z w. 


(16.b) 


Eqs. (16) constitute a measurement model and a dynamics model which are suitable for a 
simple linear Kalman filter. 

A 

Once x, y and z are estimated, the estimates of the declination angle, p , and the right 
ascension angle, a, are computed as follows 


/ A 2 A 2 

X +Y 
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/a\ 

& = tan -1 — (16.d) 

w 

Using the two-argument inverse-tangent function yields the angles in the correct quadrant 
V.2 A simple KF for (spin) phase angle determination 

It is evident from Fig. 1 that once the orientation of Z b has been determined, its 
location and the location of s in GCI coordinates define the reference for the phase angle y . 
As seen in Fig. 5, each time the sun is acquired by the sun sensor, y is equal to $ s , or $ s 
plus multiples of 2n . Therefore, the phase angle, y, has to be determined only between sun 
measurements. This is accomplished by prediction using the estimated spin rate. In order to, 
obtain superior prediction, we use a two state KF in which, during the measurement update 
stage, we improve the zero sun angle estimate and the spin rate. During the .prediction phase 
we compute the best estimate of the phase angle. The dynamics equation of this filter is 



At a sun crossing we have the following measurement 
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2n + &,=[! 



+ v„ 


(17-b) 


The value of S s is determined from the measurement s b . From Fig. 6 it is obvious that 


S s = tan 


-1 


V 




(17.c) 


z b 



Fig 6: The Geometry of 9 S 


Also at sun crossing we compute 


co„ 


2ti 


t -t, 


n-I 


(17-d) 


where t n is the present sun crossing time, and t n _, is the previous one. This "measurement" of 
the spin rate is related to the state vector in Eq. (17.a) by the measurement equation 


co. 


= [0 1 ] 


+ v„ 


(17.e) 


Consequently, at the sun crossing time, one combined measurement update is performed using 
the measurement equation 


2rt + S s 


‘1 

O' 

V 

+ 

' v r' 

- Q «" - 


0 

1_ 

CO 


_ v m _ 


(17.f) 


Once a measurement update is performed we subtract 27i from y to start a new cycle modulo 


2 7T . 
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VI. Results 

ST5 Data were used to evaluate the two filters. The data consisted of the time at which 
the measurements were taken, s b , s,, m b and m, . It should be noted that these time points 
were the instants in which the sun sensor acquired the sun. 

VL1 Spin axis filter 

We ran the spin axis filter using the following values. The value of x x , x y , and x z in . 
Eq. (16.b) was 3600 sec. The covariance matrix of the driving force noise was 
Q =diag{l0e-4 10e-4 10e-4}. The covariance matrix of the measurement noise 
corresponded to a sun sensor error of 0.1 degree, and that of the magnetometer was 10 
milliGauss, thus R = diag{3-10e-3 2.5\10e— 3}. In addition to the filtered, we also 
computed the unfiltered value of Z b in GCI coordinates. This was done using the algorithm 
presented in Eqs. (14). In Fig. 7 we present the filtered and the unfiltered components of Z b 
in GCI coordinates. The filtered values are plotted using the broken lines and the unfiltered 
values are designated by the solid lines. Obviously, the initial value of both vectors is 
identical because both were computed identically. In Fig. 8 we present the angle between the 
filtered and unfiltered unit vectors. 

It turned out that s b did not quite correspond to s, , and, similarly m b , did not quite 
correspond to m, . It was assumed that this stemmed from the fact that the measurements 
were not ideal. Another indication to this effect was the difference between the z - element of 
D^s, and of s b , and between the z - element of D^m, and of m b . This was obvious 
when checking the orthogonality of the matrix D)* computed according to Eq. (14.d). 
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Therefore we tried an additional approach to compute the filtered Z b in GCI coordinates. We 


used Eqs. (14.a - d) to compute Dj and then we applied to it Singular Value Decomposition 



Time {sj. 


Fig. 7: The Filtered and Unfiltered Components of Z b in GCI Coordinates. 



Fig. 8: The Angular Difference between the Unfiltered and Filtered Vectors. 

in order to obtain D, ort , the closest orthogonal matrix to Dj . Next we used D b l m to transform 
Sj and nij to s b and m b , respectively, and used these s b and m b as measured vectors in 
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body coordinates. The filter was then applied to the latter, and the results were termed refined- 
filtered components of Z b . It should be mentioned that the application of the TRIAD 

algorithm to the data would have also rendered an orthogonal D* ; however, it would not have 

been the orthogonal matrix closest to D* ; moreover, the result of TRIAD depends on which 

vector of the two, sun or magnetic field vector, is the one used to start the algorithm. We 
plotted the refined-filtered versus the unfiltered spin axis components in Fig. 9. The angular 
difference between the two vectors is plotted in Fig. 10. 



Time [s] 


Fig. 9: The Refined-Filtered and Unfiltered Components of Z b in GCI Coordinates. 


VI.2 Spin (phase) angle filter 

The filter described in Section V.2 was applied to the sun sensor timing data. The 
value of was 36000 sec. The covariance matrix of the white driving noise was 

Q p = diag{l0e-4 10e-4}. The initial Value of T s was computed according to Eq. (17.c), 
and the initial value of the estimated angular rate was chosen to be zero. As measurements we 
used the s b vectors which were used in the refined-filter. The covariance of the measurement 

error, r s , was chosen as 3 • lOe - 6 . The behavior of the phase angle is described in Fig. 1 1 . 
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Time Is] • 

Fig. 12: Spin Rate Estimate. * 

VII. Discussion 

As mentioned earlier, there was a discrepancy between the measured sun and the 
measured magnetic field vectors in body on one hand, and their corresponding vectors in the 
GCI coordinates on the other hand. This was obvious, when using these vectors to compute the 
transformation matrix from body to inertial coordinates or vice versa. The discrepancy 
manifested itself in non-orthogonality of the transformation Matrix. We attempted to correct 
this discrepancy by using the data to find die DCM that corresponds to each set of 
measurements, compute the orthogonal matrix closest to it, and then use it to transform the 
inertial data to body data, and treat it as the measured data to which we applied the spin axis 

filter. It is impossible to tell which plots better describe the orientation of Z b , because we do 
not have the reference attitude. 

We investigated the sensitivity of the filters. It was found that decreasing the time 
constants in the dynamics model increased the difference between the filtered and unfiltered 
components of Z b . That difference was insensitive to the change in the value of Q of that 
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model. Note that the value of R was not a design variable but was rather determined by the 
accuracy of the sun sensor and the magnetometer. The increase of Q p had very little influence 

on the estimates. Like Rof the spin axis filter, r s also was not a design variable, but was 
rather determined by the accuracy of the sun sensor. 

VHI. Conclusions 

In this work we presented a new recursive filter to estimate the attitude of a spinning 
SC. It is based on the separation of the attitude representation into the representation of the 
spin axis orientation by its components in the GCI system, and the spin (phase) angle about 
this axis. This approach enables the separation of the filter into two low-order simple filters, 
one of which estimated the slowly varying components of the spin axis, and the other 
estimated the phase angle and the spin rate. Even though the spin angle changes fast, its filter 
is very simple because the spin axis and the spin rate are almost constant. Both filters are 
independent of the SC dynamics model, which is one of the factors that make the filters so 
simple. The ambiguity problem is solved using vector calculations, which avoids tedious 
spherical geometry computations. ST5 satellite data were used to test the filters, and the 
sensitivity to filter parameter change was examined. Although there were no data to compare 
the results with, the results indicated that the filters performed well. 
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